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Abstract 

We calculate the magnetostatic energy of synthetic ferrimagnet (SyF) elements, con- 
sisting of two thin ferromagnetic layers coupled antiferromagnetically, e.g. through 
RKKY coupling. Uniform magnetization is assumed in each layer. Exact formulas as 
well as approximate yet accurate ones are provided. These may be used to evaluate 
various quantities of SyF such as shape-induced coercivity and thermal stability, like 
demagnetizing coefficients are used in single elements. 



Synthetic antiferromagnets (SAF, resp. ferrimagnets, SyF)[l, 2] consist of two thin 
ferromagnetic films of moments of same (resp. different) magnitude, strongly coupled 
antiferromagnetically thanks to the RKKY interaction through an ultrathin spacer layer, 
typically Ru 0.6 — 0.9 nm thick[3]. Hereon we consider only the case of in-plane magne- 
tized layers. SyFs are widely used to provide spin-polarized layers displaying an overall 
weak moment. One benefit is to minimize cross-talk of neighboring (e.g. memory bits) or 
stacked (e.g. in a spin-valve) elements through stray-field coupling[l, 2], such as in Mag- 
netic Random Access Memory (MRAM)[4]. SyFs are also used to decrease the Zeeman 
coupling with external fields, e.g. to increase coercivity in reference layers[5], decrease 
effects of the Oersted field in magneto-resistive or spin-torque oscillator pillars, or more 
recently boost the current-induced domain- wall propagation speed in nanostripes[6; 7]. 
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In practice SyFs are used as elements of finite lateral size. It has been shown[8] and 
it is widely used [9, 10] that for flat and magnetically soft nanomagnets of lateral size 
smaller than a few hundreds of nanometers, the macrospin approximation (uniform mag- 
netization) is largely correct. In this framework the coercive field equals the anisotropy 
field 2K/ fioMs and the energy barrier KV {V is the volume of the dot) preventing spon- 
taneous magnetization reversal equals the magnitude of anisotropy of the total magnetic 
energy £, to which all the physics therefore boils down. Elongated dots are often used to 
induce or contribute to an easy axis of magnetization and an energy barrier AS , based on 
dipolar energy. Dipolar energy is a quadratic form and thus it is fully determined by its 
value along the two main in-plane axes. For single elements A£^ = K^VAN with AN 
the difference between the two in-plane demagnetizing coefficients, and A'd = {1/2)iiqM^. 

Analytical formulas have been known for a long time to evaluate the mutual energy 
of an arbitrary set of prisnis[ll]. However while simple expressions for N and thus AS^ 
have been described, displayed and discussed for single-layer flat elements[ll. 12]. the 
analytical expressions and the evaluation of in SyFs have not been discussed in detail 
so far. Instead the studies requiring estimation of the dipolar energy in SyF. mainly 
pertaining to MRAM cells[9, 10, 13], have in the best case made use of an effective so- 
called attenuation coefficient with respect to self-energy] 10], which requires a numerical 
evaluation] 13]. The meaning and scaling laws of this attenuation coefficient have never 
been discussed in detail, hiding the physics at play. As thermal stability, coercivity and 
toggle switching fields]9, 10, 14] depend crucially on the interlayer magnetostatic cou- 
pling, it is desirable to have a simple yet accurate analytical expression for interlayer 
dipolar fields. In this manuscript we report exact analytical expressions for the magne- 
tostatics of SyFs uniformly-magnetized in each sub-layer. From the numerical evaluation 
we discuss the physics at play, while from the analytical formulas we propose an approx- 
imate yet accurate scaling law for their straightforward table-top evaluation. 

We first consider SyF prisms and name Fl and F2 the two ferromagnetic layers (Fig- 
ure 1). with magnetization aligned along z. This covers the case of both finite-size 
prisms as well as infinitely-long stripes with a rectangular cross-section. We apply for- 
mulas expressing the interaction between two parallel charged surfaces]ll], and adopt 
the convenient notation of Fijk functions, the i, j and fc-fold indefinite integrals along x, 
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Figure 1: Geometry and notations of a prismatic SyF element comprising two ferromagnetic layers Fi 
and F2. 

y and z of the Green's function i^ooo = l/r[15]. The only such function needed here is 

F220 = 7^[x{v - w)Lx + y{u - w)Ly] - xyP^ + ^r{3w - r^) (1) 

with u = , V ~ y^ , w = , r ~ y/u + v + w, = (l/2)hi[(r + x)/{r — x)] etc, 
Px ~ X eLTCtSLii{yz / xr) etc, and ~ and Px = for a; = etc. 

The integrated magnetostatic energy of a single prismatic element of thickness t is: 

fd = — V {~l)'-+''+''F22o{a5a,t6t,cS,) (2) 
<5<i,'5f ,<5ce{oa} 

which normalized to K^i yields the demagnetizing coefficient Nz. It can be verified that 
Eq. (2) coincides with the explicit formula already known[12]. The magnetostatic energy 
of a prismatic SyF element may be calculated using the same formalism . may be written 
as: 

£d = Kd,iNz{a,ti,c)Vi+KdaNz{a,t2,c)V2 

+2^IUsKd^2Nn,{a,ti,s,t2,c)^V^ (3) 

withiV„,(a,ii,s,i2,c) = ^;^^i^E5i,5,A,5ce{o,i}(-l)''^''^'"^'^ ^^220Ha,S + il<5l + 
t2<52, C(5c) is a mutual magnetostatic coefficient with a negative value, and Vi ~ atiC (resp. 
Kd,i) is the volume (resp. dipolar constant) of each single prism i. This equation of dipo- 
lar energy is a quadratic form of Mi and A/2, generalizing the definition of demagnetizing 
coefficients. 
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Figure 2; Magnetostatic energy of a SyF with c = 2a = 100 nm, Mi^2 = 10® A/m, s = 0.7 nm. (a) Sum 
of the energies of two prisms without mutual interaction, and when embedded in the SyF geometry, ti 
is kept constant at 2.5 nm, while t2 is varied, (b) Energy of the general SyF. The curved lines are those 
of minimum energy for either constant ti or t2 . The thick horizontal dotted line highlights the path for 
the SyF curve shown in (a). 

Figure 2(a) shows £d upon building a SyF via the progressive thickness increase of F2 
above Fl, considering or not the interaction between the two layers. In the latter case the 
energy increase nearly scales with , which is understandable because it is a self-energy 
(in F2 alone) . In the coupled case (a SyF) £d retains like for the uncoupled case an overall 
close-to-parabolic convex shape as can be verified with fitting, however with an initial 
negative slope. This can be understood as for low t2 the extra edge charges induced by 
an infinitesimal increase St2 mainly feel the stabilizing stray field arising from Fl, while 
for large t2 they feel more the nearby charges induced by F2 itself. Notice that, contrary 
to what could be a first guess, the minimum of £(1(^2) occurs before the compensation of 
moment (<i = i2)- This stems from the same argument as above, which is that magnetic 
charges at an edge of F2 are closer to another than to the charges on the nearby edge of 
Fl, thus for an identical amount 6t2 contribute more to £d- 

Figure 2(b) shows the full plot of £diti,t2) for s = 0.7nm. The above arguments 
appear general. From this figure let us outline three take-away messages. 1. For a given 
ti the minimum of of a SyF is found for ^2 ^ ti/2. 2. At this minimum is reduced 
by only w 20—30 % with respect to a single- layer element of thickness ti considered alone. 
3. £d roughly regains the value of the single layer at the moment compensation point 



This sheds light on results previously noticed empirically, however whose origin and 
generality had not been highlighted. Wiese et al. reported that the effective dipolar 
field anisotropy of a SyF basically scales with the inverse net moment[l6\, i.e. like the 
inverse strength of Zeeman energy. This suggests that AE^ is essentially independent of 
the imbalance of moment, which goes against the widespread belief that magnetostatic 
energy nearly vanishes upon moment compensation. Our results clarify and quantify 
this: the dipolar energy does not differ more than 20-30 % from that of a single layer 
for ^2 ^ ti (Figure 2a. dotted line). Saito et al. also reported that the thermal stability 
of Co9oFeio[3] /Ru[0.95] /Co9oFcio[5] is similar to that of Co9oFeio[3]. As explained in 
the introduction, we recall that thermal stability is determined by the energy barrier 
along the hard axis direction, with respect to the easy axis direction. In the case of 
anisotropy arising from dipolar energy and an elongated shape of the element, this bar- 
rier can be evaluated straightforwardly by calculating once along the short edge of 
the dot, and second along the long edge of the dot. Doing this we explain the findings of 
Saito et al.., whereas a reduction of 50 % would be expected on the basis of compensated 
moments (the numbers in brackets are thicknesses in nanometers). Our calculations may 
also be applicable to the cross-over of vortex versus single domain in flat disks[17] or 
vortex versus transverse domain walls in stripes[18], whose scaling law t x a = Cte may 
be derived qualitatively by equaling the energy of a vortex ^ t and that of a single- 
domain ^ a'^t(t/a) (here V ~ c?t is the volume, and tja the demagnetizing coefficient). 
Interestingly Tezuka et al. noticed that there is an optimum ferromagnetic Elm thick- 
ness at which SyAF can obtain a single-domain structure. This minimum (related to 
a minimum of demagnetization energy) is found for an imbalanced thickness in good 
quantitative agreement with Figure 2b. 

With a view to promote the use of accurate magnetostatics for SyF while eliminating 
the need for numerical evaluation, we derived approximate yet highly accurate expressions 
for £d- Figure 3a shows that to a very good approximation, is proportional to the 
width of the element (along x) and is independent of its length (along z). This is already 
accurate for a single-layer (t2 = 0), and is very accurate close to the compensation 
ti = t2 because edges then behave as lines of dipoles. whose stray field quickly decays 
with distance (~ l/?*^)- Thus boils down to a single line integral along its edge: 
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Figure 3: (a) Energy of a single layer (full symbols) and SyF (open symbols) as a function of dot length, 
i.e. along z, while a = 100 nm. (b) Energy of a single layer as a function of width (along x, open 
symbols), and length (along z, full symbols, same curve as in a), while the other in-plane dimension is 
kept constant at 100 nm. The lines are linear fits. For both plots the parameters are: Ali = lO^A/m, 
ii = i2 = 2.5 nm, s = 0.7 nm. 



Eq. (4) is the general expression, expressed in the following two lines in cartesian and 
polar coordinates (Figure 4a). E\ is the density of magnetostatic energy per unit length 
of edge, a concept once discussed in the case of single layers[19]. Equations (4-6) apply to 
an arbitrary shape of perimeter (not simply rectangles for prisms) by considering the in- 
plane angle ip between magnetization and the normal to the edge. It can be verified that 
for a SAF we have, with an accuracy better than 10 % for geometrical parameters relevant 
for practical cases, i.e. ti.2 in the range of 2 — 10 nm and sintherangeof 0.5 — 1 nm: 



The meaning of Eq. (7) is straightforward: due to the short range of interaction 




(4) 



(5) 



Ex « {l/2)Kdf 



(7) 
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Figure 4: (a) notations for the calculation of edge energy (b) integrated magnetostatic energy for various 
shapes. E is the elliptical integral of the second kind. See text for the definition of Ex. 

between dipolar lines, the density of dipolar energy is non-zero only in the vicinity of 
the edges, with a lateral range t. Thus a volume is concerned with a line density 
of energy of the order of K^i. Expressions for non-compensated cases (including single 
elements) may also be evaluated. This provides us with analytical expressions for the 
magnetostatics of SyFs for the most usual shapes (Figure 4b) . 

A scaling law sometimes used as a first guess is based on the point dipole approxima- 
tion. In this framework the energy gained by coupling Fl and F2 would roughly scale 
with K^a'^/t, resulting from two point moments MsO^t interacting like (for s t, 
and assuming lateral dimensions of the order of a). The scaling arising from our exact 
calculation is aE\ ^ K^at^ (Eq. (7) and Figure 4a). The point dipole approximation is 
thus clearly incorrect with an extra scaling (a/t)^ (see Figure 4a) which largely overesti- 
mates the dipolar coupling. This is a general argument for any flat element, where dipolar 
fields are short-ranged[20] and thus the point-dipole approximation is clearly incorrect. 

To conclude we derived exact formulas for the magnetostatics of prism SyF, and sim- 
ple yet accurate forms for SyFs of arbitrary shapes. These simple forms may be used 

straightforwardly to derive scaling laws for all aspects of SyF physics pertaining with 
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dipolar energy such as thermal stability, coercivity and anisotropy field. Notice that sim- 
ilar to the case of single flat elements edge roughness is liable to reduce significantly dipo- 
lar energy[21, 22, 23], so that the theoretical predictions need to be considered as upper 
bounds to the experimental values. The non-uniformity of magnetization is not expected 
to have a significant impact for lateral sizes below a few hundreds of nanometers [8]. 
We acknowledge useful discussions with Y. Henry, IPCMS-Strasbourg. 
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